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ABSTRACT 

A giant planet embedded in a protoplanetary disk forms a gap. An analytic 
relationship among the gap depth, planet mass Mp, disk aspect ratio hp, and 
viscosity a has been found recently, and the gap depth can be written in terms 
of a single parameter K = {Mp/M^Yh~^a~^. We discuss how observed gap 
features can be used to constrain the disk and/or planet parameters based on 
the analytic formula for the gap depth. The constraint on the disk aspect ratio is 
critical in determining the planet mass so the combination of the observations of 
the temperature and the image can provide a constraint on the planet mass. We 
apply the formula for the gap depth to observations of HL Tan and HD 169142. 

In the case of HL Tan, we propose that a planet with 0.3 is responsible for 
the observed gap at 30 AU from the central star based on the estimate that the 
gap depth is ^ 1/3. In the case of HD 169142, the planet mass that causes the 
gap structure recently found by VLA is 0.4Mj. We also argue that the spiral 
structure, if observed, can be used to estimate the lower limit of the disk aspect 
ratio and the planet mass. 
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Introduction 


High resolution imaging observations of protoplanetary disks have revealed the pres¬ 
ence of complex morphological structures in disks such as dips (e.g., Hashimoto et ah 2011), 


gap (e.g., Debes et ah 2013[ ), or spirals (e.g., Muto et ah 2012 Grady et al. 2013) 


m near 

infrared (NIR) scattered light. At submillimeter wavelengths, ALMA has revealed the ex¬ 


istence of non-axisymmetric structures in its early science operations (e.g., Casassus et al. 


2013 

Fukagawa et al.||2013 

van der Marel||2013 

Perez et al.||2014 

) and now multiple ring-like 

structures are discovered in 2014 Long Baseline Campaign ( 

ALMA Partnership 

2015 

). It 


is now possible to obtain high spatial resolution data with the beam size less than O'.'Od of 
multiple wavelengths. 

One possible origin of the detail structures in protoplanetary disks is the disk-planet 
interaction. A giant planet creates a gap and spirals in the disk along its orbit, which are 
observable with high angular resolution observations. It is therefore important to understand 
what we can learn about the physical properties of the disk and planet once we can observe 
such structures. 

In this letter, we present an estimation method of the planet mass from the gap depth 
which is potentially an observable quantity. The gap depth is related to the planet mass, the 


viscosity and the disk scale height (or disk temperature) (Duffell & MacFadyen 2013 Fung 


et al. 2014; Kanagawa et al. 2015, hereafter, K15). Using this relation, we can constrain the 


mass of a planet, if exists, by the combination of the disk image and observations of disk 
temperature and viscosity. The disk temperature can be obtained from the observation of 
optically thick gas/dust emissions. Alternatively, the spiral structure associated with the 
gap can be used to constrain the lower limit of the disk scale height. If the gap depth and 
the spiral structure are both observed, it enables us to derive the lower limit of the planet 


mass. 


In Section we present an analytic formula that describes the gap depth in terms of 
the physical parameters of the disk and the planet. In Section we apply the gap depth 
formula to the observations of HL Tau by ALMA ( ALMA Partnership 2015) and HD 169142 
by VLA (Osorio et al. 2014). In Section]^ we discuss the constraints on the disk aspect 
ratio derived from spiral structure. Section is for summary and discussion. 


2. Gap model 

In this section, we describe an analytic formula that relates the gap depth, planet mass 
and disk properties. A giant planet embedded in a protoplanetary disk creates a gap (e.g.. 
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Lin & Papaloizou 1979 

Goldreich & Tremaine 

1980 

, which has been investigated by many 

authors (e.g., Artymowicz & Lubow||1994 

Kley||1999 

Varniere et al.||2004 

Crida et al.||2006). 

Recent hydrodynamic simulations by 

Duffell & MacFadyen 

(2013 

) and 

Fung et al. 

(2014) 


give an empirical formula for the relationship between planet mass and the gap depth, which 
is derived analytically by considering the angular momentum exchange processes between 


the disk and the planet (Fung et ah 2014, and K15). We briefly explain a more accurate 
derivation given by K15, below. 

We assume a steady, axisymmetric, non-self-gravitating, geometrically thin disk. In a 
viscous accretion disk, the radial angular momentum flux Fj is given by the combination of 
the fluxes carried by advection and viscosity (|Lynden-Bell fc Pringle 1974), 


dVL 

Fj{R) = R^^Fm - 27ri?VS —, 


( 1 ) 


where i?, fl, Fm, v-, and S denote the radial distance from the central star, the angular 
velocity of the disk, the mass flux, the kinetic viscosity, and the surface density, respectively. 


Adopting the a-prescription of kinematic viscosity (Shakura & Sunyaev 1973), we write 
z/ = ah'^R^fl, where h = c/{Rfl) denotes the disk aspect ratio and c is the sound speed, 
which are dependent on the disk temperature. 

Because of a strong gravitational torque exerted by a planet on the disk, a gap opens 
along the planet orbit. We assume a sufficiently wide gap so that the torque is exerted 
mainly within a flat gap bottom with a constant surface density Sp. Then, using the WKB 
torque formula with a cutoff, the one-sided Lindblad torque, Tlb, which is exerted by the 


planet on the outer disk, is approximately given by (see Lin & Papaloizou 1979, and K15) 


Tt.r = 


2 


dR 




M, 


= 0.127r ( ^ 1 h-^RtnlEp, 


p p P’ 


( 2 ) 


where Mp and M* denotes the masses of the planet and the central star, respectively, and 
the subscript p denotes the values at the planet orbital radius Rp. In the above, we set the 
torque cutoff length A = l.ShpRp to match the one-sided torque with the results by detailed 
linear analyses of disk-planet interaction (|Takeuchi fc Miyama 1998; Tanaka et al 


2002 


Muto & Inutsuka 2009). The above assumption that the planetary torque is exerted mainly 


within the gap can be valid when the gap width is much larger than the cutoff length. 


In steady state, the mass flux is constant and the angular momentum flux radially 
changes due to the planetary torque. Thus, the angular momentum flux at the outside of 
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the gap, Re{> Rp) is equal to the sum of the angular momentum flux at Rp and the one-sided 
torque Tlb, 

Fj(i?e) = Fj(i?p)+rLB. (3) 

At the outside of the gap, the surface density is the unperturbed value Sq. The gap width 
is much smaller than Rp though the width is longer than the cutoff length. Then we can 
set Re — Rp, he — hp, and ffg — and therefor RlQeFu — R^^pFu- Under these 
approximations, substitution of Equations 0 and ([^ into (|^ yields 

SnahlR^^Ql^o = ST^ahlR^^QlEp + 0.127r (4) 

From Equation (|^, we obtain the gap depth Sp/So, as 


So 1 + 0.04A:’ 


where 


K 



( 6 ) 


(see K15 for detail derivation). Although we adopt a simple expression for the one-sided 
torque. Equation (|^ gives a good agreement with the results of hydrodynamic simulations 
as shown in Figure [Tj 


Equation ^ would be valid for K ^ 10"^. If K is appreciably large, the gap even 
becomes eccentric and non-steady, which breaks down our simple estimate (Kley & Dirksen 


2006). In the eccentric state, the gap is shallower than the prediction by Equation ([^ (|Fung 


et al. 2014) 


As seen in Equation ([^, it is related to the planet mass, the viscosity and the aspect 
ratio. Solving Equation ([^ for Mp/M*, we obtain 




5 X 10"^ 




( 7 ) 


Figure [^displays the planet mass given by Equation ([^ as a function of the gap depth Sp/So 
and the disk aspect ratio for several cases of a. The dependence of the planet mass on hp 
is strong, though a-dependence is relatively weak. The planet mass varies from 10-5 to 
~ 10 “^ if hp changes from 0.01 to 0.2 for Sp/So = 1/10 and a = lO-^. Therefore, the 
constraint on hp, or equivalently disk temperature, is important in estimating the planet 
mass. In current observations, the constraints of the disk temperature can be given, for 
example, by emissions from optically thick gas (or sometimes dust) as discussed in Section]^ 
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As discussed in Section]^ we also show that the spiral structure can constrain the lower limit 
of the disk aspect ratio. 


In application of our model to an observed gap, we assume that the distribution of dust 
is similar to the gas. Because of the dust filtration accompanied with the gap, however, the 
dust distribution can be significantly deviated from the gas , as reported by many theoretical 
studies (e.g., Paardekooper & Mellema 2004 Rice et ah 2006 Ward 2009 Zhu et ah| 2012). 
The millimeter-sized or larger dust particles strongly suffer the dust hltration. They are 
piled up at the outer edge of the gap and have a much deeper gap than the gas. Smaller 
particles with radii of ^ 0.1mm, on the other hand, are not much affected by the dust 
hltration and their distributions are similar to the gas for a gap created by a Jupiter mass 
planet (see Figure 3 of Zhu et ah 2012). Hence the above assumption in our model is valid 
for the case of such small dust particles. Also note that the critical dust size of hltration 
depends on the gas density and the turbulent viscosity. Zhu et al. (2012) did not take into 
account an additional diffusion of dust particles by the turbulence due to the Rossby wave 
instability and the Rayleigh unstable region (e.g., Zhu et al. 2014 Kanagawa et al. 2015), 
which occur at a deep gap. This additional diffusion can further increase the critical dust 
size for hltration. In application to an observed gap, it is also important to judge whether 
the hltration effect strongly alters the dust distribution at the gap or not. Such a judgement 
can be possible, by checking the dust pilling-up at the outer edge of the gap (see Section]^. 


3. Application to Observations 

In the actual observations of protoplanetary disks to date, there are several cases where 
the gap-like structures are observed. In this section, we consider two examples: HL Tau and 
HD 169142. 

observed during the 2014 Long Baseline Campaign of ALMA ( |ALMA 
. The observations show multiple rings but we focus on the ring residing 
at ~ 30 AU, which is clearly seen and reasonably resolved. 

We use the FITS image of the HL Tau that is publicly available from the ALMA 
archive site. We use the data of dust continuum emission in Band 6 and 7 in our analyses, 
and assume the distance of 140 pc. Figure |^(a) and (b) show the radial profiles of the 
brightness temperature in Band 6 and 7 between 20 AU - 50 AU averaged over the position 
angles of 130° < PA < 140°, which is close to the major axis of the disk. We use the 
position of the central brightest pixel as the location of the disk center. Since the brightness 
temperatures in Band 6 and 7 are similar at i? ^ 25 AU, the dust emission is likely to be 


HL Tau was 


Partnership 2015 
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optically thick. Therefore, we use the ratio of the brightness temperature in the two bands 
to estimate the actual gap depth. The brightness temperatures in Band 6 and 7, Thre and 
Tbrj, are given by 

Tbre = T (1 - e-"«) ( 8 ) 

and 

Tbr7 = T (1 - e-^) (9) 

where T is the disk temperature, Tq is the optical depth in Band 6 and is that in Band 7. 
The optical depths of the Band 6 and 7 are related by , where z/g = 233 GHz 

is the central frequency at Band 6 , z /7 = 344 GHz is that at Band 7, and (3 is the opacity 
index. We can derive rg and T by using Equations ([^ and ([^ if we assume jS. Figure |^(c) 
shows the radial prohles of rg in the cases of /3 = 1 and (3 = 2. By measuring the contrast of 
the maximum of r at 72 ~ 25 AU and the minimum at 7? ~ 30 AU, we estimate the depth 
of the gap to be Sp/So ~ Tmax/T-min ~ 1/3 - 1/7. 

Figure |^(d) also shows the profile of T in the cases of /3 = 1 and 2. The temperature in 
the gap region strongly depends on the assumption of /3, but that outside the gap is relatively 
insensitive to it. If we use T = 55 77 at 72 = 25 AU, and assume that the sound speed is 
c = 1 km/s(T/300 and that the central star mass is 1 Mq, we derive h ~ 0.07 as an 
estimate of the disk aspect ratio. Adopting Equation Q, we derive 


ATp 


3 - 5 X 10"^ 




( 10 ) 


if this gap is created by a planet. Note that the error in the estimate of the stellar mass 
also increases the error in the planet-star mass-ratio through hp oc c/{RpVtp) oc In 

the above estimate, it is assumed that dust particles are gas have simular distributions. 
This assumption is vaild when dust particles are too samll to be much affected by the dust 
filtration, as mentioned in the end of Section]^ We can check the validity of the assumption. 


using the observed gap structure. According to the result by Zhu et al. ( 2012 ), the dust 


surface density at the outer edge of the gap is much larger than that at the inner edge when 
the dust hltration is effective. For the gap at 30 AU of HL Tan, on the other hand, the 
distribution of the dust optical depth (see Figure [^[c]) does not show any enhance at the 
outer edge. Hence, we can judge that dust hltration is not so effective in this case. 


The second example is the disk around HD 169142. Recent observations by Osorio et al. 


(2014) of 7 mm emission shows that there may exist a gap at ~ 50 AU from the central star. 


The observations also found a knot at 50 AU, which could be regarded as compact source 
with the mass of ~ 0.6Mj. It is difficult to estimate the gap depth from their observation 


since the beam size is relatively large, but the radial prohle given in the Figure 2 of Osorio 
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et al. (2014) suggests that the gap may be as shallow as Sp/So — 0.5. We consider this 


as the upper limit of the gap depth since the model with complete gap (i.e., Sp/So = 0) 
seems to be consistent with observations. If we adopt the disk aspect ratio of ~ 0.1 and 
a = 5 X 10“^, as indicated by Osorio et al. (2014), we obtain the lower limit of the planet 
mass as ~ 0.4 Mj from Equation ([^. 


4. The Constraint on Disk Temperatnre from Spiral Structures 


We have seen that the constraint on disk temperature, or disk aspect ratio, is important 
in inferring the planet mass from the gap structure. As discussed above, the disk temperature 
is obtained from the disk emissions. Here, we discuss an alternative constraint on disk 
aspect ratio coming from a spiral density wave, which is a natural consequence of disk- 
planet interaction. The spiral structures have been observed in infrared observations so far 


e.g., Muto et al. 2012, Grady et al. 2013 Currie et al. 2014). We expect such features may 


be observed at sub-millimeter wavelengths in near future when better sensitivity at higher 
spatial resolution is available. 

This shape of the spiral density wave induced by disk-planet interaction is related to 


the disk aspect ratio (or temperature) and is given in an analytic form (Rahkov 2002 Muto 


et al.||2012 ). The shape of spiral waves excited by a planet is determined by four parameters, 
Rc,0q,6 and hp (see eq. [1] Muto et al. 2012). The launching location of spiral waves (very 


close to the planet location) is denoted as {Rc,0o)- The parameter 6 denotes the profile of 
sound speed c oc R~^. The disk aspect ratio h changes as h = hp {R/RpY^^^^"^. 

For demonstration, we have performed two-dimensional hydrodynamic calculations us¬ 


ing FARGO (Masset 2000) to obtain the gap and spiral structures caused by a Jupiter mass 
planet at Rp = 40 AU in a disk with a = 10“^, hp = 0.1, and 5 = 0.25. Figure |^(a) and 
(b) are the optically thin dust thermal emission at 230GHz obtained by the simulation {k 
is assumed to be 2 x 10“^ cm^/g per unit gas mass and the dust and gas are assumed to be 
well-mixed), and the radial cuts of the emission along the wave and in the direction opposite 
to the planet. Figure |^(c) and (d) are the same as (a) and (b), but the emission is convolved 
with a circular Gaussian beam with an FWHM of 0.065” (the distance to the disk is set as 
150 pc). Even in the convoluted image, the emission of wave is about 10% larger than the 
emission from non-wave region near the planet. Such spiral feature would be detectable with 
future ALMA or TMT observations. We ht the spiral structure, and the resulting y^-map 
(Figure]^ indicates that we can constrain the lower limit of hp to be ~ 0.07. The gap depth 
can be estimated by Fp/So ~ 1/5 from Figure]^ Therefore, from Equation (j^, we derive 
the lower limit of the planet mass as Mp/M* > 4x 10 “^(q!/ 10“^)^/^. This lower limit actually 

































coincides with the actual planet mass within a factor of ~ 2.5 in this particular case. 


5. Summary and Caveats 

In this letter, we have shown an analytic derivation of the depth of the gap created by a 
planet embedded in a protoplanetary disk. Equation Q gives the relationship between the 
gap depth, planet mass, disk temperature (aspect ratio) and viscosity. Since the gap depth 
strongly depends on the disk aspect ratio, the constraints on the disk temperature is critical 
in estimating the planet mass from the gap structure. The disk temperature is obtained by 
the emissions from optically thick gas/dust, as shown in Section]^ We have also argued 
that the spiral wave can constraint the disk temperature in the observation of optically 
thin emission. Our model assumed that dust particles and gas have similar distributions. 
This assumption is vaild when dust particles are too small to be much affected by the dust 
filtration. The assumption would be verified if the dust piling-up is not observed at the 
outer edge of the gap. 

We have applied Equation ([^ to the dust continuum observations of HL Tau and 
HD 169142 and derived the estimates of the mass of the planet, which could be respon¬ 
sible for the observed gaps (at 30 AU of HL Tau and at 50 AU of HD 169142). At the gap of 
HL Tau, the dust filtration does not seem so effective becasue no significant dust pilling-up 
is seen at outer edge of the gap. 

We also note that in both the observations of HL Tau and HD 169142, the gaps we 
have focused on are only marginally resolved. Therefore, the gap depth quoted in this paper 
should be regarded as upper limit. Future high resolution observations by ALMA and/or 
TMT should put good constraints on the gap depth. When complex morphologies are found 
in protoplanetary disks, simple analytic formulae should be useful for interpreting (at least 
a part of) such structures and derive physical quantities. 


KDK thanks Tatsuya Takekoshi for helpful comments. This paper makes use of the 
following ALMA data: ADS/JAO.ALMA^2011.0.00015.SV. ALMA is a partnership of 
ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRG 
(Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The 
Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work is partially 
supported by JSPS KAKENHI Grant Numbers 23103004, 26103701, 26800106. Numerical 
computations were carried out on Cray XC30 at Center for Computational Astrophysics, 
National Astronomical Observatory of Japan. This work is also greatly benefited from 
Workshop hosted by Institute of Low Temperature Science, Hokkaido University ’’Recent 


Development in Studies of Protoplanetary Disks with ALMA”. 
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Fig. 1.— The gap depth as a function of the parameter K. The solid line denotes Equa¬ 
tion (|^. The crosses, circles and triangles represent the gep depth obtained by hydrody¬ 


namic simulations done by Duffell & MacFadyen (2013), Fung et ah (2014) and Varniere 


et ah (2004), respectively, (entire runs of iF < 8000 in these papers are shown) 
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Fig. 2.— Planet mass as a function of the gap depth and the disk aspect ratio. The 
relationship between these values are given in Equation Q and three panels correspond to 
the cases with a = 10“^ (left), 10“^ (center) and 10“^ (right). 
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Fig. 3.— (a): Observed radial profile of the brightness temperature of dust continuum 
emission in HL Tan disk along the major axis. The data is averaged over PA=135° ± 5°. 
Solid line shows the Band 6 data and the dashed line shows the Band 7 data, (b): The ratio 
of brightness temperature in Band 6 and 7. (c): The radial prohles of the optical depth, (d): 
The temperature. The derivations of (c) and (d) are described in Section]^ 
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Fig. 4.— (a): The gap and spiral structure caused by a Jupiter mass planet in the disk 
with a = 10“^, hp = 0.1 and 6 = 0.25. The black dashed line indicates the best-£t curve 
for the spiral shape and the white dashed line indicates the spiral used for htting. (b): The 
radial cuts of emission along the spiral (dashed) and at a direction opposite to the planet 
(solid), (c) and (d): The same as (a) and (b), but the emission is convolved with a circular 
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Fig. 5.— The of the htting to the spiral structure. 
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